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ABSTRACT 

The extraordinary properties of nonlinear optical propagation processes in double-domain positive/negative 
index metamaterials are reviewed. These processes include second harmonic generation, three- and four-wave 
frequency mixing, and optical parametric amplification. Striking contrasts with the properties of the counterparts 
in ordinary materials are shown. We also discuss the possibilities for compensating strong losses inherent to 
plasmonic metamaterials, which present a major obstacle in numerous exciting applications, and the possibilities 
for creation of unique ultracompact photonic devices such as data processing chips and nonlinear-optical sensors. 
Finally, we propose similar extraordinary three-wave mixing processes in crystals based on optical phonons with 
negative dispersion. 

Keywords: Negative-index plasmonic metamaterials, backward electromagnetic waves, three- and four-wave 
mixing processes, coherent energy exchange, optical parametric amplification, compensating losses, optical 
phonons to replace negative optical magnetism, photonic microdevices. 

1. INTRODUCTION 

In the late 1960s, V. G. Veselago considered the propagation of electromagnetic waves in an fictitious, isotropic 
medium with simultaneously negative dielectric permittivity e and magnetic permeability /i and showed that it 
would exhibit very unusual properties!^! Specifically, the simultaneously negative dielectric permittivity, e < 0, 
and magnetic permeability, \x < 0, would lead to a negative refraction index and to a left-handed triplet of 
the electric field, magnetic field and the wavevector. The energy flow (Poynting vector) in this case appears 
counter-directed with respect to the wavevector. This is rather counter-intuitive and in a sharp contrast with 
normal, positive-index materials (PIMs), also referred to as right-handed materials (RHMs). Negative- index 
materials (NIMs) , also referred to as left-handed materials (LHMs) , do not exist naturally. It was also generally 
accepted in physical optics that the magnetization at optical frequencies is negligible and, hence, did not play 
any essential roleP In accordance with this, the magnetic permeability \i was normally set to be equal to one in 
the basic Maxwell's equations describing linear and nonlinear optical processes. 

Metamaterials, i.e., artificially designed and engineered materials, can have properties unattainable in nature, 
including a negative refractive index. As outlined above, the backwardness of electromagnetic waves, i.e., the 
phenomenon of counter-directed energy flow and phase velocity for electromagnetic waves, does not exist in natu- 
rally occurring materials. Its realization has become feasible owing to the advent of nanotechnologies. Significant 
progress in the design of bulk metamaterials has been achieved during the last few years (see, e.g.J^I). This 
dictates a revision of many concepts concerning electromagnetic propagation processes in man-made materials, 
which hold the promise of revolutionary breakthroughs in microwave and photonic device technologies. 7 

Negative refraction can be implemented to develop a wide variety of devices with enhanced and uncommon 
functions, such as hyper-resolution lenses and optical cloaking devices. Nanostructured metamaterials are ex- 
pected to play a key role in the development of all-optical data processing chips. The opposite directions of 
the energy flow and wavevector (phase velocity) in NIMs determines unique linear and nonlinear optical pro- 
cesses. Initially, metamaterials with a magnetic response and a negative refractive index were fabricated for the 
microwave and terahertz frequency ranges ISH121 optical frequency range imposes increasing difficulties and 
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challenges for metamaterials. A negative magnetic permeability in the optical range, which is a precursor for a 
negative refraction, has been demonstrated in RefsP^E^ Finally, simultaneously negative magnetic permeability 
and electric permittivity in the optical telecommunicationfrequency range, and a hence negative refractive index, 
were experimentally demonstrated for metal-dielectric metamaterials!^^ For review, see RefP Great progress 
was recently achieved in fabricating plasmonic NIMs, including chiral NIMs, and is described in RefePH^^SD 

To satisfy the casuality principle, a NIM must be lossy PsED The majority of NIMs realized to date contain 
metal nanostructures. These structures introduce strong losses inherent to metals that are difficult to avoid, 
especially in the visible range of frequencies. Irrespective of their origin, losses constitute a major hurdle to 
the practical realization of the unique optical applications of these structures. Therefore, developing efficient 
loss-compensating techniques is of a paramount importance. So far, the most common approach to compen- 
sating losses in NIMs is associated with the possibility of embedding amplifying centers in the host matrix. 
The amplification is supposed to be provided through a population inversion between the energy levels of the 
embedded centers. Such a method to remove absorption losses and dissipation by optical amplification, e.g. by 
introducing a lasing medium instead of a passive dielectric into a negative-positive index composite and to utilize 
alternating amplifying layers instead of a bulk NIM, was suggested inPSl Although there are problems with this 
approach (for example, regarding amplified spontaneous noise, the influence of surface plasmons and the extreme 
localization of fields in optical amplification), it appears that it could furnish improved resultsPESfHD R am an 
amplification can be also used for removing the loss obstacleP-3 Exciting progress was achieved during the last 
year in producing population-inversion-based amplification in plasmonic metamaterialspSEU including recent a 
breakthrough with compensating losses in NIMsP^ 

The main emphasis in the studies of NIMs has been placed so far on linear optical effects. Nonlinear optics 
in NIMs, especially frequency mixing, still remains a less-developed branch of photonics. It has been shown that 
NIMs, which include structural elements with non-symmetric current-voltage characteristics - its "meta-atoms," 
can possess a nonlinear magnetic response at high frequencies PUSH Such meta-atoms have highly controllable 
magnetic and electric responses. Hence, NIMs can combine unprecedented linear and nonlinear electromagnetic 
properties. The feasibility of crafting NIMs with nonlinear optical (NLO) responses in the optical wavelength 
range, such as second-harmonic generation (SHG) and third-harmonic generation in magnetic metamaterials, 
has been experimentally demonstrated pBESl Qn a fundamental level, the NLO response of nanostructured 
metamaterials is not yet completely understood or characterized. Nevertheless, it is well established that local- 
field-enhanced nonlinearities can be attributed to plasmonic nanostructures. Most negative-index (NI) plasmonic 
metamaterials are metal-dielectric composites made of metallic nanostructures shaped as split rings, horseshoes, 
fishnets, or chiral structures that are plunged into a dielectric host. This makes ah initio calculations of nonlinear 
propagation processes a complex task. Such works fall into a special, rapidly expanding research category. 
Recent progress in studying NLO responses of plasmonic metamaterials, both in theory and proof-of-principle 
experiments, is described in RefsPsHMl anc j references therein. 

The important properties of SHG in NIMs in the constant-pump approximation were discussed in RefP^ for 
loss-less, semi-infinite materials and in RefsJSKD f or a s i a b f a finite thickness. Unlike ordinary NLO materials, 
the latter appeared to bring major qualitative effects. The possibility of exact phase-matching for waves with 
counter-propagating energy-flows has been shown in RefP^l for the case when the fundamental wave falls in the 
negative-index frequency domain and the second-harmonic (SH) wave lies in the positive-index domain. The 
possibility of the existence of multistable nonlinear effects in SHG was also predicted in RefP^l As outlined here, 
absorption is one of the most challenging problems that needs to be addressed for the realization of practical 
applications of NIMs. A transfer of the near-field image into aSH frequency domain, where absorption may be 
lower, was proposed in REfsP^EH as a possible means to overcome dissipative losses and thus enable the superlens. 
The propagation of microwave radiation in nonlinear transmission lines, which are the one-dimensional analog 
of NIMs, was investigated in RefP^ The possibility to achieve conversion efficiency up to 12% was shown in 
RefP^ for phase-matched SHG in a strongly absorbing, bulk NIM. In that case, the input picosecond-pulse peak 
intensity was assumed to be about 500 MW/cm 2 , the nonlinear susceptibility of the metamaterial was assumed 
to be x 2 ~ 80 pm/V, and the attenuation depth of the NIM was about 50 /mi. The possibility of achieving phase- 
and group- velocity matching was shown for the above-indicated pulses based on the Drudc dispersion model for 
metals. Nonlinear-optical processes in NIMs possess unusual, sometimes counter-intuitive, properties. They are 
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described below with the example of SHG in loss-less NIMs of a finite thickness taken from our recent paperPHl 
The extraordinary NLO properties of such propagation processes in NIMs, including three- wave mixing (TWM), 
four- wave mixing (FWM) and optical parametric amplification (OPA), also have been predicted and are in stark 
contrast with their counterparts in natural materials EsHU] Striking changes in the properties of nonlinear pulse 
propagation and temporal solitonsj^l spatial solitons in systems with bistabilityj^HM] g a p solitonsj^ and optical 
bistability in layered structures including NIM o 68 ^ 69 ! have been revealed. A review of some of the corresponding 
theoretical approaches is given in Refs.^E^ Frequency-degenerate multi-wave mixing and self-oscillations of 
counter-propagating waves in ordinary materials have been extensively studied earlier because of their easily 
achievable phase matching. Phase matching for TWM and four-wave mixing (FWM) of contra-propagating 
waves that are far from degeneracy seem impossible in ordinary materials and presents a technical challenge in 
the metamaterials. It has become achievable only recently due to the advances in nanotechnology.^H^ The 
possibility and extraordinary properties of TWM, including mirrorless self-oscillations, was proposed in RefP^ 
(and references therein) more than 40 years ago (see also Refs.^^) for two co-propagating waves with nearly- 
degenerate frequencies that fall within an anomalous dispersion frequency domain and gives rise to generation 
of a far-infrared difference-frequency counter-propagating wave in an anisotropic crystal. However, far-infrared 
radiation is typically strongly absorbed in crystals, which presents an unavoidable and strong detrimental factor. 
The proposal was not realized until recently. For the first time, a TWM backward-wave (BW) mirrorless optical 
parametrical oscillator (BWMOPO) with all three significantly different optical wavelengths was realized in 
Reff^ Phase- matching of counter-propagating waves has been achieved in a periodically poled NLO crystal with 
the period on a submicrometer scale. Both in the proposal and in the experiment J^l the opposite orientation of 
wavevectors was required to establish a distributed feedback and to produce mirrorless OPO. This was due to the 
fact that ordinary, positive-index crystals were proposed for the implementation. As outlined, a major technical 
problem in creating BWMOPO stems from the requirement of phase matching for the opposite orientation of 
wavevectors in PIMs. This paper reviews the counter-intuitive effects and unusual features in the energy exchange 
between coupled ordinary and contra-propagating backward electromagnetic waves associated with three types of 
nonlinear optical processes: second harmonic generation beyond the constant fundamental wave approximation, 
three-wave mixing, and four-wave mixing. 

The extraordinary properties of second-harmonic generation in a frequency double-domain positive/negative 
matamaterial are studied in Section [21 Both semi-infinite and finite-length NIM slabs are considered and com- 
pared with each other and with ordinary PIMs for SHG. The feasibility of a nonlinear-optical mirror converting 
the incident radiation into a reflected SH is shown and specific features of the process are investigated both for 
continuous- wave and pulsed regimes. The paper is organized as follows. The "backward" features of the electro- 
magnetic waves in NIMs are discussed in Section [27T1 The basic equations for the negative-index fundamental 
and positive-index SH waves are derived and the relevant Manley-Rowe relations are discussed and compared 
with those in PIMs in Section [272] The unusual spatial distributions of the field intensities for SHG in a NIM slab 
of a finite thickness are described in Section 12.41 for the model of a loss-free, double-domain metamaterial and 
a continuous-wave regime. These distributions are compared with their counterparts in ordinary PIMs (Section 
12 . 3[) and semi-infinite NIM slabs ( Section 12.51) . Section [2751 also analyzes the properties of a nonlinear-optical 
mirror with a semi-infinite thickness that converts an incident fundamental beam into a reflected SH one. The 
role of absorption and phase mismatch is discussed in Section 12.61 The extraordinary features associated with 
SHG in the pulsed regime, which include significant differences in the pulse shapes for the fundamental and SH 
radiations dependent on the intensity and on the pulse duration of the incident pulses, are considered in Section 
I2~T1 

The extraordinary features of coherent, NLO energy transfer from the control optical field(s) to the coupled 
contra-propagating negative-phase (negative index, NI) and positive-phase (positive index, PI) waves through 
three- and four- wave mixing and related optical parametric amplification are described in Sections [3] and [U The 
uncommon phenomenon of generating a contra-propagating wave at an appreciably different difference frequency 
in the direction of reflection is investigated. The feasibility and unparalleled features of such energy transfer are 
shown, which stem from the unusual fact that the energy flow of one of the coupled electromagnetic waves is 
contra- directed relative to the others, whereas the wavevectors for all coupled waves remain parallel. Such an 
opportunity makes phase matching of counter-propagating waves much easier and is offered by the backwardness 
of electromagnetic waves that is natural to NIMs. Consequently, distributed-feedback features become possible, 
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while the antiparallel orientation of the wavevectors of the coupled waves is no longer required. The different 
properties attributed to such an all-optically tailored nonlinear optical mirror-amplifier for the case of a negative- 
index control field and for the alternative option of a negative-index idler are discussed and compared in Sections 
13.11 and 13.21 The characteristic magnitudes of the required intensity of the control field and of the slab thickness 
are given in Section 13.31 

The feasibility of independently crafting of a negative index and a resonantly enhanced four-wave mixing 
nonlinearity through embedded four-level centers is considered in Section 0] The characteristic values of the 
required density of the embedded resonant nonlinear-optical centers, their relaxation properties, and the intensity 
of the control field are given in Section |4~T1 The possibility of substituting a nanostructured negative-index metal- 
dielectric composites, which requires sophisticated fabrication techniques, with extensively used and studied 
ordinary crystals in order to simulate the unparalleled properties of coherent NLO energy exchange between 
the ordinary and backward waves is described in Section [5] Finally, the possibilities of implementing originally 
strongly absorbing microscopic samples of plasmonic, nanostructured, metal-dielectric composites for remote, 
all-optical tailoring of the transparency and reflectivity of metamaterial films as well as the concepts of several 
options of unique, ultracompact photonic devices based on the outlined three- and four-wave processes are 
summarized in the concluding Section [6j 

2. SECOND-HARMONIC GENERATION IN A LOSSY, DISPERSIVE 

DOUBLE-DOMAIN SLAB 

2.1 Wavevectors and Poynting vectors for the fundamental and second-harmonic waves 
in a double-domain metamaterial 

We consider a loss-free material that is left-handed at the fundamental frequency uj\ (ei < 0, //i < 0), whereas it 
is right-handed at the SH frequency uj 2 — 2wi (e 2 > 0, /x 2 > 0). The relations between the vectors of the electric, 
E, and magnetic, H, field components and the wavevector k for an electromagnetic wave written as 

E(r,i) = E (r)exp[-i(urt-k-r)] + c.c., (1) 
H(r,i) = Ho(r)exp[-i(wt-k-r)] + c.c., (2) 

traveling in a loss-free medium with dielectric permittivity e and magnetic permeability \x are given by the 
equations 

k x E = (w/c)/iH, k x H = -(cj/c)eE, (3) 

ViE(r,t) =-^H(r,t), (4) 

which follow from Maxwell's equations. Equations ([3]) show that the vector triplet E, H and k forms a right- 
handed system for the SH wave and a left-handed system for the fundamental beam. Simultaneously negative 
material parameters (e < and [i < 0) result in a negative refractive index n = —y/JIi. As seen from Eqs. (fTj) 
and ([2|), the phase velocity v p ^ is co-directed with k and is given by v p h = (k/fc)(w//c) = (k/fc)(c/|n|), where 
k 2 = n 2 {w/c) 2 . In contrast, the direction of the energy flow (Poynting vector) S with respect to k depends on 
the signs of e and /i: 

S(r,t) = -^[E x H] = -r^— [H x k x H] - ^H 2 = -^E 2 . (5) 

v ' 47r 47rwe 47rwe inufj, 

Since we have taken the material to be loss-free, we assume here that all indices of e, fj, and n are real numbers. 
Thus, the energy- flow Si at uji is directed opposite to ki, whereas S 2 is co-directed with k 2 . 

2.2 Basic equations and the Manley-Rowe relations for SHG process in NIMs and PIMs 
2.2.1 Basic equations 

We assume that an incident flow of fundamental radiation Si at uj\ propagates along the z-axis, which is 
normal to the surface of a metamaterial. According to ([5|), the phase of the wave at u>i travels in the reverse 
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Figure 1. The difference in the phase-matching geometry and in the intensity distribution for the fundamental and the 
second-harmonic waves, /i? and h%, between a slab of the left-handed (negative- index) metamaterial, (a), and a right- 
handed material (ordinary, positive- index) material, (b). 



direction inside the NIM [Figfjja)]. Because of the phase- matching requirement, the generated SH radiation 
also travels backward with energy flow in the same backward direction. This is in contrast with the standard 
coupling geometry in a PIM [Figfjjb)]. As shown below, the basic features of the process associated with electric 
or magnetic nonlinearities are similar. Following Ref.p^ we assume that a nonlinear response is primarily 
associated with the magnetic component of the waves. Then the equations for the coupled fields inside a NIM 
in the approximation of slow-varying amplitudes acquire the form: 

dA 2 /dz = ia 2 Alexp(-iAkz), (6) 
dAx/dz = ia 1 A 2 A* 1 exp(iAfcz). (7) 

Here, Afc = k 2 — 2ki\ o\ — (eiw 2 /fcic 2 )87rx' 2 \ o 2 = (e2wf/fc2C 2 )47rx^; x is the effective nonlinear sus- 
ceptibility; A 2 and A\ are the slowly varying amplitudes of the waves with the phases traveling against the 
z-axis, 

Hj(z,t) — Aj exp[— i(tOjt + kjz)] + c.c, (8) 

j = {1, 2}; uj 2 = 2u\; and fci. 2 > are the moduli of the wavevectors directed against the z-axis. We note that 
according to Eq. (Q}, the corresponding equations for the electric components can be written in a similar form, 
with ej substituted by /ij and vice versa. The factors [ij were usually assumed to be equal to one in similar 
equations for PIMs. However, this assumption does not hold for the case of NIMs, and this fact dramatically 
changes many conventional electromagnetic relations. 
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2.2.2 The Manley-Rowe relations 

The Manley-Rowe relationship for the field intensities and for the energy flows follow from Eqs. ([5]) - ([7]): 



t\ dz 2t2 dz dz dz 

The latter equation accounts for the difference in the signs of t\ and £2, which brings radical changes to the 
spatial dependence of the field intensities, as discussed below. 

In order to outline the basic difference between the SHG process in NIMs and PIMs, we assume in our 
further consideration that the phase-matching condition fc 2 = 2fcj is fulfilled. The spatially-invariant form of the 
Manley-Rowe relations follows from equation ©: 

|^i| 2 Ai + l^2| 2 A2 = c ! (10) 

where C is an integration constant. With t\ = —62 Equation (fTOf takes the form: 

\A 1 \ 2 -\A 2 \ 2 = C, (11) 



Equations (|10[) and (JTTJ) predict a concurrent decrease of the amplitudes of both waves along the z-axis so that 
the difference between the squared amplitudes remains constant through the sample, as schematically depicted 
in Fig. Qla). Hence, Eqs. (TTOj) and (jlip also predict that the difference between the number of pairs of photons 
fiwi and the number of photons fvjj2 remains constant through the sample. This is in striking contrast with 
the requirement that the sum of the squared amplitudes and the sum of the corresponding photon numbers is 
constant in the analogous case in a PIM, as schematically shown in Fig. Hfb) and described in textbooks on 
nonlinear optics. 

Ultimately, the described extraordinary features of SHG in a double-domain NIM and their dependence on 
the slab geometry stem from the backwardness of electromagnetic waves in NIMs and the opposite directions of 
the energy flows in fundamental-harmonic (FH) and SH waves. As shown in the next sections for three-wave 
mixing processes, different depletion rates caused by different absorption indices at frequencies of the contra- 
propagating wave may qualitatively change their propagation properties, the distribution across the slab, and 
the output characteristics depending on the specific NLO propagation process. 

We now introduce the real phases and amplitudes of the fields as Ai$ = exp(i0i,2)- Then the equations 
for the real amplitudes and phases 

dh-^/dz — U2h\sin^ , 
dhi/dz = — o~ihih2sin 1 fy , 

d^/dz = -{2aih2-(T2h\/h2)cos^, (12) 
* = 2 (z)-2<Mz)., 

which follow from Eqs. © and (UJ, show that if any of the fields becomes zero at any point, the integral (|10p 
corresponds to the solution with the constant phase difference 2(j)\ — §2 — tt/2 over the entire sample. 

2.3 SHG in PIMs 

The equations for the slowly varying amplitudes corresponding to the ordinary coupling scheme in a PIM, shown 
in Fig. [TJb), are readily obtained from Eqs. (JSJ) - © by changing the signs of ki and k 2 . This does not change 
the integral (jlOl) ; more importantly, the relation between e\ and £2 required by phase matching now changes to 
£1 = £2, where both constants are positive. The phase difference remains the same. Because of the boundary 
conditions ft-i(O) = h w and h 2 (0) — h 2 o — 0, the integration constant becomes C — h\ . Thus, the equations for 
the real amplitudes in the case of a PIM acquire the form: 

h 1 (z) = ^hi -h 2 (zy, (13) 

dh2/dz = K[hl ~h 2 (z) 2 }, (14) 
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with the known solution 



h 2 (z) = /i 10 tanh(z/z ), (15) 
h\(z) = hi /cosh(z/z ), z = [nh w }~ 1 . (16) 



Here, k = (i-2^2 1 '^2 c2 )47rXg*' / ■ The solution has the same form for an arbitrary slab thickness with decreasing 
fundamental and increasing SH squared amplitudes along the z-axis, as shown schematically in Fig. [TJb). 

2.4 SHG in a NIM slab 

Now consider phase-matched SHG in a loss-less NIM slab of a finite length L. Equations ([6]) and (fTTj) take the 
form: 

hi(z) 2 = C + h 2 (z) 2 , (17) 
dhi/dz = -k[C + h 2 (z) 2 }. (18) 

Taking into account the different boundary conditions in a NIM as compared to a PIM, hi(0) = h\o and h^{L) = 0, 
the solution to these equations is 



h 2 = VCtaxi[vCK{L-z)], (19) 
hi = Vc/cob[VCk(L- z)], (20) 

where the integration parameter C depends on the slab thickness L and on the amplitude of the incident funda- 
mental radiation as 

VCkL = cos" 1 (VC/hio). (21) 

Thus, the spatially invariant field intensity difference between the fundamental and SH waves in NIMs depends on 
the slab thickness, which is in strict contrast with the case in PIMs. As seen from Equation (|17p . the integration 
parameter C = h\{z) 2 — h2(z) 2 now represents the deviation of the conversion efficiency 77 = ^20/^10 from unity: 
(C/hf ) = 1 — 7]. Figure H] shows the dependence of this parameter on the conversion length zq = (k/ii ) -1 . The 




Figure 2. The normalized integration constant C /h\ and the energy conversion efficiency n vs. the normalized length of 
a NIM slab. 



figure shows that for a conversion length of 2.5, the NIM slab, which acts as nonlinear mirror, provides about 
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80% conversion of the fundamental beam into a reflected SH wave. Figure [3] depicts the field distribution along 
the slab length. One can see from the figure, with an increase in slab length (or intensity of the fundamental 
wave), the gap between the two plots decreases while the conversion efficiency increases (comparing the main 
plot and the inset). 




Figure 3. The squared amplitudes for the fundamental wave (dashed line) and SHG (solid line) in a loss-less NIM slab of 
a finite length. Inset: the slab has a length equal to one conversion length. Main plot: the slab has a length equal to five 
conversion lengths. The dash-dot lines show the energy conversion for a semi-infinite NIM. 



2.5 SHG in a semi-infinite NIM 

Now we consider the case of a semi-infinite NIM at z > 0. Since both waves disappear at z — > oo due to the 
entire conversion of the fundamental beam into SH, C = 0. Then equations p7|) and (fT5|) for the amplitudes 
take the simple form 

h 2 (z) = hx{z), (22) 
dh 2 /dz = -nh\. (23) 

Equation (|22|) indicates 100% conversion of the incident fundamental wave into the reflected second harmonic at 
z = in a loss-less, semi-infinite medium provided that the phase-matching condition Afc = is fulfilled. The 
integration of ([23| with the boundary condition hi(0) = hio yields 

(z/z ) + 1 

Equation (|24p describes a concurrent decrease of both waves of equal amplitudes along the z-axis; this is shown 
by the dash-dot plots in Fig. |3l For z > zo, the dependence is inversely proportional to z. These spatial 
dependencies, shown in Fig. [3l are also in striking contrast with those for the conventional process of SHG in a 
PIM, which are shown in various textbooks [compare, for example, with Figfljb)]. 
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2.6 SHG in a lossy dispersive NIM slab 

It is convenient to introduce effective amplitudes, Oj, and nonlinear coupling parameters, Xj, which are defined 

as 

a J = \[WJkj A h X = \f%k2/\l4»2\4*Xc8- 

Here, Xog is the effective nonlinear susceptibility, and the quantities \a,j | 2 are proportional to the photon numbers 
in the energy fluxes. 

In the general case of SHG in a dispersive double-domain material with a(wi) ^ a(uj2) and Ak — fc 2 — 2fci ^ 0, 
the equations for the amplitudes ejj take the form: 

dai/dz = — i1Xa\a2 exp(iAfcz) — (a\/2)a\, (25) 
da 2 /dz = iXa\ cxp(-iAfcz) + (a 2 /2)a 2 . (26) 

Here, ai j2 are the absorption indices at the corresponding frequencies. The equations account for opposite 
signs of ei, [i\ and e 2 , /i 2 . Three differences distinguish these equations from their counterparts in ordinary 
PIMs - opposite signs with the NLO coupling factors, opposite signs with the absorption indices, and boundary 
conditions to be applied to the opposite sides of the slab. These differences stem from the backwardness of the 
fundamental wave in the NIM in this case. Note that, for a NIM with an electric nonlinearity, the corresponding 
equations for the electric components can be written in a similar form, with ej substituted by fij and vice versa. 




Figure 4. Effect of absorption (a) and phase mismatch (b) on the distribution of the FH and SH fields across the slab and 
on the output values. 

Equations (l25l) - (|26l) do not have analytical solutions. Figures [4] and [5] depict the results of a numerical 
analysis of the chief features of SHG in a dispersive, absorptive slab for boundary conditions a\{z = 0) = aw, 
a 2 (z = L) = 0. The quantity T\ — |ai(z)| 2 /|aio| 2 represents the distribution of the fundamental beam inside the 
slab determined by both absorption and energy conversion as well as the ultimate slab transparency, T\{z = L). 
The quantity 7/ 2 = 2|a 2 (z)| 2 /|aio| 2 represents the quantum conversion efficiency of the process, g = Xaio, and 
g^ 1 is the nonlinear conversion length, i.e., the characteristic slab thickness that is required for significant energy 
conversion and is dependent on the input intensity of the fundamental field. The plot for u\ — a 2 = in Fig. [4fa) 
displays two curves with a constant gap, which determines the output value of about 80% for the reflected beam 
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a L=2.3, a L=1, Akl_=7t a L=2.3, a L=1 , Akl_=7t 




Figure 5. Dependence of the distribution of the FH and SH fields across the slab and of their output values on the intensity 
of the input fundamental beam for AkL — n. 



(z = 0) at 2wi for the given input intensity of the fundamental field at u)\. By comparing this plot with the 
others, we see that absorption causes variable a gap that increases in the output area z > 0, so that the output 
value of SHG drops to about 40% for a±L — 2.3 and a 2 L — 1 and the same intensity of the input fundamental 
beam. For a\L — 2.3 and <zio — > 0, the transparency of the slab at u>i reaches up to 10%, which is characteristic 
for plasmonic samples. Figure EJb) shows that the conversion efficiency decreases with an increase in phase 
mismatch, however it is less than it would be for the ordinary coupling schemes in PIMs. FigurelUa) shows that, 
for the given parameters, the SH field concentrates in the area close to the exit facet of the slab at z = 0. The 
dependence on g at z — indicates that a significant part of the incident beam can be converted into SH for 
gL 1 even for the given phase mismatch and absorption rates. This is shown in Fig. [5jb), which depicts the 
significant depletion of the FH field with an increase of its input intensity and the slab length (parameter gL). 

2.7 Pulsed SHG in a transparent NIM slab 

The equations for the quantities a\ and a 2 introduced above are 

—i2ga*a 2 exp (iAkz) 2~ ai 

iga\ exp (— iAkz) + ~~^ a 2- (27) 

Here, |<3, lj2 1 2 are the above-introduced slowly varying amplitudes that are proportional to the photon numbers 
in the energy fluxes corresponding to the pulse maximum, and Vi are the group velocities for the corresponding 
pulses. After introducing £ = z/L and r = t/ At, d — L/v\At, where d is the slub thickness L reduced by the 
input pulse length and At is the duration of the input fundamental pulse, Eqs. (1271) take the form: 

—i2ga\ai exp (iAfc£) 2~ ai 

igaj exp (-iAfc£) + y ai (28) 



1 dai da\ 

vi dt dz 

1 da,2 <9q,2 

vo dt dz 



da\ dai 



dr 
, da 2 
dT 



da 2 
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The input pulse shape is taken as being close to a rectangular form 

F(t) = 0.5 ( tanh T ° + * - T - tanh ^f^] , (29) 
\ or St J 

where St is the duration of the pulse front and tail, and r is the shift of the front relative to t = 0. The 
magnitudes St — 0.01 and tq — 0.5 have been selected for numerical simulations. 
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(d) 



Figure 6. Input T\(z = 0) and output T\(z = L) pulses of fundamental and negative- index SH radiation 7/2(2 = 0) for 
relatively short pulses of fundamental radiation, (a) Input pulse area S10 = 0.9750; output pulse areas S\l = 0.5031, 
S20 = 0.2392. (c) Input pulse area S10 = 0.9750; output pulse areas Sil = 0.0396, S20 = 0.4742. 

The unusual properties of SHG in NIMs in the pulsed regime stem from the fact that it occurs only inside the 
pulse of fundamental radiation. Generation begins on the leading edge of the pulse, grows towards its trailing 
edge, and then exits the pulse with no further changes. Since the fundamental pulse propagates across the 



11 



gl_=2, d=1 , a 1 =a 2 =0, Ak=0 



gL=2, d=1 , a =a =0, AkL=2n 








(a) 

gl_=1 0, d=1 , a, =a n =0, Ak=0 

» 1 2 







--T1 2 (Z=0) 






--^(2=0) 

—2^(2=1.) 








1 ' 
1 ' 
1 f 

1 f 


1 

\ 

\ 

\ 






1 

0.8 
0.6 
0.4 
0.2 

0. 







(b) 

gL=10, d=1, a =a =0, AkL=2n 









--2t! 2 (z=0) 




1 ,1 
' >■ 

1 1 




--^(2=0) 




1 : 1 
I 

' 1 


: 


—21/2=1) 






■ 

/ V 

\ V 

—/ r 





(c) 



(d) 



Figure 7. Input 7i(,z = 0) and output T\(z = L) pulses of fundamental and negative- index SH radiation 7/2(2 = 0) for 
relatively long pulses of fundamental radiation, (a) Input pulse area S10 = 0.9900; output pulse areas Sil = 0.2516, 
S20 = 0.3692. (c) Input pulse area S10 = 0.9900; output pulse areas Sil = 0.0161, S20 = 0.4870. 
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Figure 8. Effect of phase mismatch on SHG in a N1M slab for the case of strong, relatively long pulses. 
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Figure 9. Effect of difference in the group velocities on SHG in a NIM slab for the case of strong, relatively long pulses. 
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slab, the duration of the SH pulse is longer than the fundamental one. Depletion of the fundamental radiation 
along its pulse length and the conversion efficiency depend on its initial maximum intensity and phase matching. 
Ultimately, the overall properties of SHG, the pulse length, and the photon conversion efficiency depend on the 
ratio of the fundamental pulse and slab length. This extraordinary behavior is illustrated in Figs. (SJIHl 

Figures HHS] display the input shape of the fundamental pulse Ti = |oi(z)| 2 /|aio| 2 at z = when its leading 
front enters the medium and the results of numerical simulations for the output fundamental pulse when its 
tail reaches the slab boundary at z = L as well as the shape and the conversion efficiency of the output second 
harmonic pulse rj 2 = |a 2 (z)| 2 /|aio| 2 traveling against the z-axis, when its tail passes the slab's edge at z — 0. For 
clarity, here, the medium is assumed loss-less and group velocities of the fundamental and SH pulses assumed 
equal for Figs. 151151 

Figure [6] corresponds to the fundamental pulse four time shorter than the slab thickness. It shows increase 
of the conversion efficiency with increase the intensity of the input pulse. It is followed by the shortening of the 
SH pulse. Phase mismatch causes changes in the depletion and shape of the output fundamental pulse. However 
overall conversion rate does not change significantly. The outlined properties satisfy to the conservation law: the 
number of annihilated pair of photons of fundamental radiation (Sio — «Sil)/2 is equal to the number of output 
SH photons S^o- 

Figure [7| displays the corresponding changes for a longer input pulse length equal to 0.5 of the slab thickness. 
Here, conversion efficiency increases at a lower input intensity because of the longer conversion length. The 
changes in the SH pulse length and conversion efficiency with increasing input intensity appear less significant. 

Figure [8] shows the effect of phase mismatch on the shape of the SH pulse for relatively long pulses and a 
high-intensity input pulse of fundamental radiation. 

Figure [9] shows that the properties of the output SH pulse do not change significantly with an increase of the 
group velocity difference up to 20% for long pulses and a high-intensity fundamental field. 

3. THREE-WAVE MIXING OF CONTRA-PROPAGATING ELECTROMAGNETIC 
WAVES: TAILORED TRANSPARENCY AND REFLECTIVITY OF A PLASMONIC 
METASLAB, AND APPLICATIONS TO ULTRACOMPACT NONLINEAR-OPTICAL 

DATA PROCESSING CHIPS AND SENSORS 

An extraordinary electromagnetic property of NIMs stems from the fact that the energy flow and phase velocity 
of electromagnetic waves become counter-directed inside the NIM slab. Usually, a negative refractive index exists 
only inside a certain frequency band. The metamaterial remains ordinary, PI, outside that band. Consider a slab 
of thickness L that possesses a quadratic nonlinearity, which enables three- wave mixing processes W3 = u>i + 10-2 ■ 
The amount of coherent energy transfer between the coupled light waves increases with the decrease of phase 
mismatch AkL, where Ak = fc 3 — k 2 — fci . This dictates the requirement of co-directed wavevectors for all coupled 
waves. Consequently, the energy flux for the wave that falls in the negative-index frequency domain appears 
contra-directed relative two others. Figure [TU] depicts three corresponding possible options for the phase-matched 
TWM NLO coupling of ordinary and backward waves. 

For the continuous- wave regime, the equations for the slowly varying amplitudes of the coupled waves can be 
written in the form 

dai/dz = ±[iX±a2a3 exp(iAkz) — (a?i/2)ai], (30) 
daijdz = iX2Cila3 exp(iAkz) — (a?2/2)a2, (31) 
da^/dz = ^[1X30102 exp(— iA/cz) — (a3/2)a3]. (32) 

The minus and plus signs refer to the negative- and positive-index waves, correspondingly. Here, Ak = ks—k2—ki, 
and all three waves experience strong dissipation described by absorption indices 0^1,2,3 ■ The slowly varying 
effective amplitudes of the waves, Q>e,m.ji (j — {1?2,3}) and nonlinear coupling parameters, Qe,mi 

for the electric 
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(a) (b) (c) 

Figure 10. Three different phase matching options for three-wave mixing of ordinary and backward light waves, (a) Si, 2 
and ki,2 are energy fluxes and wavevectors for the ordinary, positive index signal and generated idler; S3 and k3 - for the 
negative index control field, (b), (c) Alternative schemes, (b) The NLO chip amplifies the negative- index signal traveling 
against the control beam [n(oji) < 0] and frequency upconverts it to the contra-propagating beam, (c) The NLO chip 
generates a negative-index idler and, therefore, shifts the frequency and reflects back a positive-index signal traveling 
along the control beam [n(ai2) < 0]. 



(e) and magnetic (m) types of quadratic nonlinearity can be conveniently introduced as 

a ej = yj\ej/kj\Ej,X e = ^\k 1 k 2 /€ 1 € 2 \2TrXe j ) ; 
a m j = \J\nj/kj\H jy X m = v^iV^i/t^TTx^-. 

The quantities |aj| 2 are proportional to the photon numbers in the energy fluxes. Equations for the amplitudes 
Oj are identical for the both types of the nonlinearities. 

We note the following three fundamental differences in Equations (|30| and (|32|) as compared with their 
counterparts in ordinary, PI materials [Eq. (pTTj)]. First, the signs with gi or g% are opposite to that with g 2 . 
This is because the corresponding ej < and < 0. Second, the opposite sign appears with a>j because the 
corresponding energy flow Sj is against the z-axis. Third, the boundary conditions for the negative-index wave 
are defined at the opposite side of the sample as compared to the ordinary waves because their energy flows are 
counter-directed. These modifications lead to dramatic changes in the equation solution as compared with the 
exponential dependence on z that is characteristic for the ordinary, naturally occurring, PI materials. 

Consider the example depicted in Fig.fTUta). We assume that the wave at with the wavevector k x directed 
along the z-axis is a PI (ni > 0) signal. Usually it experiences strong absorption caused by metal inclusions. The 
medium is assumed to possess a quadratic nonlinearity v/ 2 -* and is illuminated by the strong, higher-frequency 
control field at 0J3, which falls into the NI domain. Due to the TWM interaction, the control and signal fields 
generate a difference- frequency idler at uj 2 = u>a — u>\, which is also assumed to be a PI wave (n 2 > 0). The 
idler, in cooperation with the control field, contributes back into the wave at u>i through the same type of TWM 
interaction and thus enables OPA at ui\ by converting the energy of the control fields into the signal. In order 
to ensure effective energy conversion the induced traveling wave of nonlinear polarization in the medium and 
the coupled electromagnetic wave at the same frequency must be phase matched. Hence, all phase velocities 
(wavevectors) must be co-directed. Since n{uj^) < 0, the control field is a backward wave, i.e., its energy flow 
S3 = (c/47r)[E3 x H3] is directed against the z-axis. Such a device can be employed as NLO sensor. It can 
be conveniently and remotely interrogated to actuate frequency up-conversion and amplification of a signal 
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directed towards the remote detector. Such a signal could be incoming far-infrared thermal radiation emitted 
by the object of interest or signal that carries important spectral information about the chemical composition 
of the environment, for example. Two other schemes depicted in Fig. UOf b'l.fc) offer different advantages and 
operational properties for NLO NIM-based devices. The research challenge is that such an unprecedented NLO 
coupling scheme leads to changes in the set of coupled nonlinear propagation equations and boundary conditions 
compared to the standard ones known from published literature. This in turn results in dramatic changes in 
their solutions and in the multi-parameter dependencies of the operational properties of the proposed devices. 

Unusual propagation and energy-conversion properties in double-domain NIM/PIM slabs are readily seen, 
for example, in the coupling scheme of Fig. ITUT a) and with a loss-free medium. Then the following Manley-Rowe 
relations can be derived from Maxwell's equations for the slowly varying amplitudes: 



Here, \a,j\ present photon numbers. Equation (|33[) predicts that the difference of the numbers of photons ftwi 
and hiui remains constant through the sample, which indicates their creation in pairs due to the split of photons 
fiu>3. However, Eqs. (|34|) predict that the differences of the numbers of photons tvjj\ and Hu)3 as well as of huj2 and 
ht03 also remain constant through the sample. This looks like a breaking of the energy conservation law and is in 
seemingly striking difference with the fact that the sum of the corresponding photon numbers is constant in the 
analogous case in a PIM. Actually, such unusual dependencies stem from the fact that the waves propagate in the 
opposite direction. Consequently, extraordinary distributions of these fields across the slab and the dependence of 
their output values on the linear and nonlinear optical properties of the given NIM and on the input intensities 
of the coupled fields are expected, especially when the conversion efficiency becomes large. Particularly, the 
conversion rate is expected to grow across the slab with a different rate than in an ordinary medium with a 
standard coupling geometry. Equations (|34l) indicate unusual feedback, which provides correlated depletion of 
the control field on one hand and growth of the signal and the idler on the other hand, so that the difference must 
remain constant along the metaslab. Absorption would change this behavior, which may strongly depend on 
the absorption dispersion and on the phase mismatch. Investigation of the indicated dependencies is important 
for optimizing the operational properties of the proposed sensor. Consequently, accounting for the boundary 
conditions that must be defined at the opposite facets of the slab, this suggests extraordinary distributions of 
these fields across the slab and extraordinary dependencies of their output values on the linear and nonlinear 
optical properties of the given NIM and on the input magnitudes of the fields. Such dependencies are essentially 
different for the coupling schemes depicted in Fig. [TO] 

3.1 Properties of a nonlinear optical reflector and amplifier tailored by the 
negative- index control field 

Consider the scheme depicted in Fig. [TU^a) . An analytical solution is not possible for the problem that underlies 
the concept of the proposed sensor. Figures [TTJ and [12] display the results of numerical simulations for this 
case. Here, z is the length across the slab of thickness L, T 3 (z) — \a 3 (z) / o 3 l\ 2 and T\(z) = |ai(z)/a3i| 2 are 
transparencies, r)2(z) — \a 2 {z) / a^ 2 is the photon conversion efficiency, g — Xaa(L), X e j = yj \ ki k 2 / e i e 2 1 ^ n xfj ; 
and a$L = o,3{L) is the amplitude of the input control field. Absorption indices, otj, for the coupled field 
are indicated, and Afc = 0. Figure [TT] illustrates the case of a weak input signal so that the depletion of 
the control field due to the conversion becomes significant only in the vicinity of z=L. Figures lllf a) and (b) 
show the possibility to achieve many orders of amplification of the signal traveling against the control beam 
and its conversion to the frequency-shifted wave for the intensity of the incoming control field corresponding 
to gL ss 15. ..20. Then the numbers of the output photons fiwi and hio 2 make about 10% of that of the input 
contra-directed control field, which means amplification on the order of 10 7 . Figure [T2l demonstrates a stronger 
energy-conversion effect due to a higher input intensity of the signal which, however, leads to lower overall 
amplification. 



d(|ai| 2 -M 2 )/^ = 0, 
d(|a 3 | 2 - \ ai \ 2 )/dz = 0, d(\a 3 \ 2 ~ \a 2 \ 2 )/dz = 0. 



(33) 
(34) 
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Figure 11. Distribution of the coupled fields across the slab and output characteristics of the amplified transmitted, 
generated, and depleted control fields for the case of a weak input signal a in = 10~ 4 a,3L- 



3.2 Tunable NLO mirror employing a negative-index idler: weak signal and idler 
approximation 

The physical principles of the proposed nonlinear-optical micromirror, which also can be viewed as an optical 
data-processing chip, are based on difference-frequency generation of a backward, NI wave in a strongly absorbing 
double-domain NIM slab [Fig. [T3lah which is equivalent to Fig. HOfc)]. Simultaneously, the incident PI signal 
propagating through the slab experiences OPA. Only the reflected-wave frequency is assumed to fall in the 
negative-index frequency domain. Here, two incident co-propagating PI waves would produce an ail-optically 
controlled reflectivity accompanied by a tunable frequency shift of the reflected beam. A strong control field 
at W3 and a strong incident wave at u>2 are both assumed to be PI. The generated difference-frequency wave at 
uj\ = CJ3 — UJ2 is NI and, therefore, a backward wave. All three waves experience strong dissipation described by 
absorption indices £¥1,2,3- 

First, consider a simplified model that neglects depletion in the control field due to the TWM conversion 
and absorption. Hence, the field is assumed to be uniform across the slab. This model allows one to understand 
the basic properties of TWM coupling of backward waves and difference-frequency generation of a backward NI 
wave. For the case of a spatially homogeneous, constant control field and a real, nonlinear susceptibility, an 
analytical solution to Eqs. (|5D ]) -(|5T j) can be found, and then the reflectivity, r\ = |ai(0) / 'a^l 2 , is given by the 
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Figure 12. Distribution of the coupled fields across the slab and output characteristics of the amplified transmitted, 
generated, and depleted control fields for the case of a strong input signal aio — O.I03L. 



equation 



n 



(g/R) sin RL 



cos RL + (s/R) sin RL 



(35) 



whereas the corresponding transmission factor for the PI incident wave, T2 = |a2(£)/<J2o| , is 



T 2 = 



exp{-[(a 2 /2)-s]L} 



cos RL + (s /R) sin RL 



(36) 



Here, L is the thickness of the slab, g = go = Xa 30 , R = \J ' g 2 — s 2 , and s = [(o>i + c tg)/ ^ — [iAk/2]. It is seen 
that the reflectivity and transmittance experience a set of "geometrical" resonances*^! For example, at s = 0, 
ri=tan 2 (gi), r\ = l/cos 2 (gL) and they tend to infinity at gL — > (2j + l)n/2, (j=0, 1, 2, ...), which indicates 
the possibility of mirrorless self-oscillations. A similar behavior is characteristic for distributed-feedback lasers 
and is equivalent to a great extension of the NLO coupling length. It is known that even weak amplification 
per unit length may lead to lasing provided that the corresponding frequency coincides with high-quality cavity 
or feedback resonances. Such a giant enhancement of the conversion efficiency in the vicinity of "geometrical" 
resonances is in striking contrast with the exponential dependencies T2 cx exp(2gL) known for the counterparts 
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Figure 13. (a) Coupling scheme: S2 - incident beam, S3 - control beam, Si - generated negative-index (reflected) beam, 
(b) Characteristic dependence of reflectivity, n, and transmittance, T2, of the NLO mirror on the product of nonlinear 
susceptibility, strength of the control field and thickness of the NIM slab with account for absorption of all coupled fields. 
Here, reflectivity can be switched from zero to magnitudes exceeding 100%. 

in ordinary nonlinear materials. Figure ITST b'l displays such resonances in an initially strongly opaque slab. It 
shows the dependence of the reflectivity r± and transmittance T2 of the slab on the parameter gL for specific 
absorption indices at the frequencies of the coupled waves indicated in the figure and Ak = k% — &2 — &i = 0. 
It is seen that the modulation properties of the reflected and transmitted beams differ greatly. The transmitted 
beam experiences amplification for any magnitude of the parameter gL above a certain level. On the contrary, 
the intensity of the reflected beam at uj\ can be varied over a wide range by changing the intensity of the control 
field. It can be modulated from zero to magnified values exceeding generation threshold for both coupled waves. 
The Manley-Rowe equations for the given coupling scheme predict the sum of |ai| 2 and \d2\ 2 to be conserved 
across the slab. Note that in a transparent PI medium, the difference of these values is invariant across the 
slab. It appears that absorption of the control field causes broadening the resonances but does not destroy 
the resonance behavior. The transmission minima are dependent on the ratio of absorption rates; this is in 
contrast with the reflectivity minima, which remain robust. Ultimately, the NLO coupling under consideration 
may provide for significantly sharper growth of the entangled counter-propagating photons in the vicinity of the 
"geometrical" resonances than the conventional schemes for the same parameters of the slab and the control 
field. The described extraordinary properties can be utilized for sensing, filtering and conversion of weak light 
signals. 

Alternatively, the phase mismatch causes a decrease of the reflectivity maxima and an increase of the minima. 
Reflectivity becomes relatively robust against phase mismatch with an increase of the intensity of the control 
field. It drops dramatically in the range of small phase mismatch and then remains relatively robust at the 
lower level within the range of greater phase mismatch, as seen in Fig. 1141 The outlined properties of the NLO 
mirror are determined by the interplay of several processes that have a strong effect on the NLO coupling of 
contra-propagating waves and, consequently, on their distributions inside the slab. 

The depletion of the control field due to energy conversion to the reflected and transmitted beams may have 
a strong effect on the processes under investigation. Particularly, it limits the amplification in the maxima and 
may even change the resonance behavior. Basically, the reflected wave has a different frequency than the incident 
beam. The simulations predict that the quantum conversion efficiency with respect to the control field may be up 
to several tens of percent, which indicates great enhancement of the reflected and transmitted beams with respect 
to the incident, positive-index, weak field and the reflectivity that may significantly exceed 100%. Ultimately, 
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a U2.3, a L=1,a =0 0^1=2.3, a 2 L=1 ,a L=2.1 




Figure 14. Reflectivity vs. intensity of the control field and phase mismatch for different absorption indices for the coupled 
waves. 



the simulations show the possibility of all- optical tailoring and switching of the reflectivity and transparency of 
such a mirror over a wide range by changing the intensity of the control field. More details regarding the given 
coupling scheme and the one displayed in Fig. fTUT b) can be found in p^fglaSllaSBaHEi] 

3.3 Estimates 

The characteristic magnitude of the parameter gL, which is required to realize the effects predicted above, is on 
the order of 1. Assuming \^ ~ 10~ 6 ESU (~ 10 3 pm/V), which is on the order of that for CdGeAs2 crystals, 
and a control field of I ~ 100 kW focused on a spot of D ~ 50 /xm in diameter, one can estimate that the 
typical required value of the parameter gL ~ 1 can be achieved for a slab thickness in the microscopic range of 
L ~ 1/im, which is comparable with that of the multilayer NIM samples fabricated to date. 
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4. QUANTUM ENGINEERING OF NONLINEARITY AND COHERENT 

QUANTUM CONTROL 



Figure [TS] depicts the basic principles of independent quantum engineering of nonlinearity and its coherent 
quantum control. Here, a resonantly enhanced, higher-order, electrical x NLO response and four-wave mixing 




(a) (b) (c) (d) 

Figure 15. Coupling geometry and alternative schemes of four-wave mixing through embedded, resonant, nonlinear-optical 
centers, (a, c) Coupling geometry for four-wave mixing of backward and ordinary electromagnetic waves. Si, ki and uii 
are the energy flux, wavevector and frequency for the backward- wave signal; S2, k2 and 102 - for the ordinary idler; S3, 4, 
k3,4 and 0*3,4 - for the ordinary control fields, (b, d) Corresponding alternative schemes of quantum controlled four-wave 
mixing in embedded, resonant, nonlinear-optical centers with different ratios of the signal and idler absorption rates 
and nonlinear susceptibilities, (b) Shortest-wavelength negative-phase signal where, depending on the partial relaxation 
rates, parametric amplification can be assisted by the idler's population-inversion or Raman-type amplification, (d) 
Longer-wavelength negative-phase signal where, depending on the partial relaxation rates, parametric amplification can 
be assisted by the signal's incoherent amplification attributed to population-inversion or Raman-type gain. 

are employed in a composite metamaterial with embedded NLO centers (quantum dots, ions, or molecules). 
Two different options correspondingly depicted in panels (a), (b) and (c), (d) offer different means of quantum 
control by the driving fields at frequencies W3 and W4. Only the wave at lo\ is presumed to be a negative- index 
one, and all others are PI waves. Four-wave mixing u\ — lj^ + co^ — 0J2 replaces the three-wave mixing process 
described above, so that in this case two strong fields are used for controlling the process. Resonant processes at 
the depicted quantum transitions driven by the strong control fields make all linear and nonlinear characteristics 
of the composite intensity-dependent. Among the specific features attributed to the resonant coupling is the fact 
that the nonlinear susceptibility becomes complex in the vicinity of the introduced resonances, which is followed 
by a phase shift between the fields and the polarizations. Due to the different contributions of the populations 
of the coupled energy levels to Xi an d to x^\ the latter becomes different and intensity-dependent.^ The 
ratio of real and imaginary parts of xf^ varies with the change of the resonance offsets. It is also different for 

schemes depicted in Figs.[T57b) and (d). It appears that the detrimental role of even a large phase mismatch can 

(3") 

be diminished or eliminated through the appropriate adjustment of the real and imaginary parts of \) ■ Such 
a possibility is not available for off-resonant wave mixing. The properties of the coherent energy transfer from 
the control fields to the negative-phase signal depend strongly on the ratio of the absorption rates and NLO 
susceptibilities for the signal and the idler. Figures I15f b) and (d) present two alternative options for quantum 
control of NLO propagation processes. Figure [TSJb) depicts the scheme that offers the possibility of incoherent 
amplification (population inversion or Raman-type) for the idler ISMHUIZI This possibility depends on the set of 
quantum relaxation rates inherent to the specific doping agent. Figure 115( d) presents a scheme that offers an 
alternative possibility. Here, the idler frequency is close to a higher-frequency transition from the ground state, 
and the signal corresponds to a lower-frequency transition between the excited states. Hence, absorption for the 
idler is typically larger than for the signal. This situation appears advantageous for robust transparency tailoring. 
No incoherent amplification is possible here for the idler, whereas incoherent absorption for the signal can be 
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controlled and turned to amplification. Such a possibility is also contingent on the appropriate partial population 
and coherence relaxation rates attributed to the embedded centers. In all of the schemes outlined above, the 
linear and nonlinear local parameters can be tailored through quantum control by varying the intensities and 
frequency-resonance offsets for combinations of the two control driving fields. 

Depending on the partial transition relaxation rates associated with levels m and <?, the numerical model 
relevant to the coupling scheme depicted in Fig. I15f d) offers neither the possibility of population inversion nor 
Raman-type amplification. The fact that all involved optical transitions are absorptive determines essentially 
different features of the overall resonant loss-compensation technique compared to that proposed in RefsP^E^^ 
The following model, which is characteristic of ions and some molecules embedded in a solid host, has been 
adopted: energy level relaxation rates T n = 20, T g = T m = 120 (all in I0 6 s _1 ); partial transition probabilities 
Ign = 50, 7„ m = 90, (all in I0 6 s _1 ); homogeneous transition half-widths r s ;=1.8, r TOn =1.9, r gn =l, r m ;=1.5, 
T mg = 5 x 10~ 2 , r ni = 5 x 10~ 3 (all in I0 11 s" 1 ); Ai = 756 nm and A 2 = 480 nm. The density-matrix methooP 
has been used for calculating intensity-dependent local parameters while accounting for the quantum nonlinear 
interference effects. This allows one to investigate the changes in absorption, amplification, and refractive indices 
as well as in the magnitudes and signs of the NLO susceptibilities caused by the control fields. These changes 
depend on the population redistribution over the coupled levels, which in turn strongly depends on the ratio of the 
partial transition probabilities. Numerical simulations have led to promising conclusions regarding the enhanced 
reflectivity and amplification of the incident beam as well as possibilities of modulating these characteristics of 
the proposed microscopic metadevices.^'^^'^HsU 

4.1 Estimates 

For transition properties that are characteristic of molecules embedded in a solid environment, estimates show 
that the required intensity of the control fields at 103 and W4 are on the order of / ~ IW/(0.1mm) 2 . With a 
resonant absorption cross-section 040 ~ 10~ 16 cm 2 , which is typical for transitions with an oscillator strength 
of about one, and a concentration of embedded centers N ~ I0 19 cm~ 3 , one estimates aio ~ I0 3 cm -1 and the 
required slab thickness in the microscopic range L ~ (1 — 100)/iin. The contribution to the refractive index by 
the impurities then is estimated as An < 0.5(A/47r)a4o ~ 10~ 3 , which essentially does not change the negative 
refractive index. 

5. BACKWARD- WAVE OPTICAL PHONONS AND DISTRIBUTED 
NONLINEAR-OPTICAL FEEDBACK 

A different scheme of TWM for ordinary and backward waves has been proposed previously.^ It builds on stim- 
ulated Raman scattering (SRS) where two ordinary electromagnetic waves excite a backward, elastic, vibrational 
wave in a crystal and thus initiate TWM. The possibility of elastic BWs was predicted by L. I. Mandelstam in 
1945,^21 who also pointed out that negative refraction is a general property of BWs. The basic idea of replacing 
the negative-index composites with ordinary crystals and mimicing the uncommon properties of coherent NLO 
energy exchange between the ordinary and backward waves was shown in RefP^ and is reviewed below. 

The dispersion curve uj(k) of phonons in crystals containing more than one atom per unit cell has two branches: 
acoustic and optical. For the optical branch, the dispersion is negative in the range from zero to the boundary of 
the first Brillouin zone (Fig. [TBI . Hence, the group velocity of optical phonons, vg r , is antiparallel with respect 
to its wavevector, k£ , and phase velocity, v£ . This is because v gr — duj(k)/dk < 0. Optical vibrations can 
be excited by light waves through two-photon (Raman) scattering. This gives us the grounds to consider these 
crystals as analogues to a medium with a negative refractive index at the phonon frequency and to examine 
the processes of parametric interaction of three waves. Two of these waves are ordinary electromagnetic waves. 
The third wave is the backward wave of elastic vibrations with directions of its energy flow and wavevector that 
are opposite to each other. We investigat the extraordinary nonlinear propagation and output properties of the 
Stokes electromagnetic wave in one of two different coupling geometries, depicted in Fig. 1161 with both utilizing 
the backward elastic waves. Such unusual properties are in contrast with those attributed to the counterparts in 
the standard schemes that build on the coupling of co-propagating photons and phonons. They are also different 
from the properties of the phase-matched mixing of optical and acoustic waves for the case when the latter has 
its energy flux and wavevector directed against those of one of the optical waves. The revealed properties can 
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Figure 16. Negative dispersion of optical phonons and two phase-matching options: (a) co-propagating, and (b) contra- 
propagating fundamental, control, and Stokes signal waves. Insets: relative directions of the energy flows and the wavevec- 
tors. 



be utilized for the creation of optical switches, filters, amplifiers and cavity-free OPOs based on ordinary NLO 
crystals without the requirement of periodically poling at the nanoscale!^ Therefore, Fig. [16] depicts the concept 
of substituting a sophisticated, nanostructured, negative-index metal-dielectric composite by extensively used 
and studied ordinary crystals in order to mimic the unparalleled properties of coherent NLO energy exchange 
between ordinary and backward waves. 

6. CONCLUSION 

In conclusion, we have described the extraordinary properties of coherent, nonlinear-optical processes in negative- 
index metamaterials, such as second-harmonic generation, difference- frequency generation and optical parametric 
amplification. The feasibility of compensating the strong losses inherent to nanostructured, negative-index, 
plasmonic, metal-dielectric composites is outlined and illustrated with the aid of numerical simulations. The 
elimination of the detrimental role of optical losses in negative-index metamaterials is the key problem that 
limits the numerous revolutionary applications of this novel class of electromagnetic metamaterials. All-optical 
tailoring of transparency and reflectivity is shown to be possible through coherent, nonlinear-optical energy 
transfer from the ordinary control wave to a negative-index, backward-wave field. This property is intrinsic to 
negative-index metamaterials. It is shown that besides the nonlinearity attributed to the building blocks of the 
negative-index host, a strong nonlinear optical response in the composite can be independently provided through 
embedded, resonant, four-level nonlinear-optical centers. This opens the way to independent nanoengineering 
and adjustment of the negative index and nonlinearities of metamaterials. In addition, we have described the 
opportunity for quantum control over the local optical parameters of the metamaterial in this case, which 
employs constructive and destructive quantum interference tailored by two auxiliary driving control fields. Such 
a possibility is proven with the aid of a realistic numerical model. The possibility of mimicing the extraordinary 
properties of coherent energy transfer between ordinary and the contra-propagating, backward waves inherent 
to negative-index, plasmonic metamaterials is described by making use of an easily available class of crystals 
with negative dispersion for optical phonons. Among the possible applications of the described nonlinear-optical 
processes are a novel class of miniature, all-optically tailored sensors, mirrors, frequency-tunable narrow-band 
filters, quantum switches, amplifiers, cavity-free microscopic optical parametric oscillators that allow generation 
of entangled counter-propagating left- and right-handed photons, and all-optical data-processing chips. The 
unique, unparalleled features of the underlying processes are outlined and include the strongly resonant behavior 
with respect to the material thickness, the density of the embedded resonant centers and the intensities of the 
control fields, the feasibility of negating the linear phase-mismatch introduced by the host material, and the role 
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of absorption or, conversely, the supplementary nonparametric amplification of the idler. 
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